Lifting of nodes by disorder in extended-s state superconductors: application to 

ferropnictides 
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We show, using a simple model, how ordinary disorder can gap an extended-s [Aig) symme- 
try superconducting state with nodes. The concommitant crossover of thermodynamic properties, 
particularly the T-dependence of the superfluid density, from pure power law behavior to an acti- 
vated one is exhibited. We discuss applications of this scenario to experiments on the ferropnictide 
superconductors. 
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I. INTRODUCTION 

When a new class of unconventional superconduc- 
tors is discovered, it is traditional to try to deter- 
mine the symmetry class of the order parameter, as 
this may provide clues as to the nature of the pair- 
ing mechanism. Direct measurements of the super- 
conducting order parameter are not possible, and indi- 
rect phasc-scnsitivc measurements are sometimes diffi- 
cult because they typically involve high-quality surfaces^. 
On the other hand, a set of relatively straightforward 
experimental tests are available to determine the ex- 
istence of low-energy quasiparticle excitations and, in 
some cases, their distribution in momentum space. As 
in earlier experiences with other candidate unconven- 
tional systems like the high-Tc cupratcs and heavy 
fermion materials, the symmetry class of the newly 
discovered ferropnictide superconductors^ is in dispute 
at this writing, due in part to differing results on 
superfluid densitjs^iii^i^iiiSii^, angle-resolved photoemis- 
sion (ARPES) "'^°'^^i^^'^'^i^^'^^ , nuclear magnetic resonance 
(NMR)i^iiiii^ii^, Andreev spectroscopy20i2ii22^^ 
other probes. In some cases, results have been taken 
to indicate the absence of low-energy excitations, i.e. a 
fully developed spectral gap. In others, low-energy power 
laws have been taken as indication of the existence of or- 
der parameter nodes. It is not yet clear whether these 
differences depend on the stoichiometry or doping of the 
materials, or possibly on sample quality. 

In parallel, microscopic theoretical calculations of the 
pairing interaction in the ferropnictide materials have at- 
tempted to predict the momentum dependence of the or- 
der parameter associated with the leading superconduct- 
ing instability. Using a 5-orbital parameterization of the 
density functional theory (DFT) bandstructure, Kuroki 
et al^ performed an RPA calculation of the spin and 
orbital contributions to the interaction to construct a lin- 
earized gap equation. They found that the leading pair- 
ing instability had s-wave (Aig) symmetry, with nodes 



on the electron-like Fermi surface ( "/? sheets" ) , and noted 
that the next leading channel had dx^^y2 (-Big) symme- 
try. Wang et al^ studied the pairing problem using the 
functional renormalization group approach within a 5- 
orbital framework, also finding that the leading pairing 
instability is in the s-wave channel, and that the next 
leading channel has dx2_y2 symmetry. For their inter- 
action parameters, they found however that there were 
no nodes on the Fermi surface, but there is a significant 
variation in the magnitude of the gap. Other approaches 
have obtained Aig gaps which change sign between the 
hole and electron Fermi surface sheets but remain ap- 
proximately isotropic on each shee t^^'^^ . 

We also recently presented calculations of the spin and 
charge fluctuation pairing interaction within a 5-orbital 
RPA framework^^, using the DFT bandstructure of Cao 
et ali^ as a starting point. Our results indicated that 
the leading pairing channels were indeed of s (Aig) and 
dx2-y2 symmetry, and that one or the other could be 
the leading eigenvalue, depending on details of interac- 
tion parameters. We also gave arguments as to why these 
channels were so nearly degenerate, and pointed out some 
signiflcant differences in the states compared to those 
found by Kuroki et ali^. Finally, we noted that, within 
our treatment and for the interaction parameter space 
explored, nodes were found in all states, generally on the 
a sheets for the c?-wave case and the /? sheets for the s- 
wave case, but that in the latter case the excursion of the 
order parameter of sign opposite to the average sign on 
the sheet was small, and might be lifted by disorder—. 

It is the purpose of this paper to explore the pos- 
sibility that the nodes of an extended-s state of the 
type discussed in Ref. [2^ are lifted by disorder, and 
consider the consequences for experimental observables. 
In the interest of simplicity, we initially neglect many 
complicating aspects of the problem, in particular the 
multisheet nature of the Fermi surface, and focus pri- 
marily on the sheet found in each case (s or d) which 
has nodes. In this case the problem is similar to one 
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which was studied earher in the context of potential 
extended-s states in a single-band situation for cuprate 
superconductors^2i^i22. For the pnictides. current inter- 
est centers around the isotropic sign-changing extended-s 
state proposed by Mazin et al^, one where a gap is mo- 
mentum independent over two independent Fermi surface 
sheets, but has a different sign on each. In this case, it is 
known that only intcrband scattering is pairbrcaking, be- 
cause it mixes the two signs and suppresses the overall or- 
der parameteri22i^i2ii^. This effect has been claimed to 
account for low-energy excitations observed in NMR^i^^ 
and superfluid density^ experiments. In general, impu- 
rities with screened Coulomb potential will have a larger 
intraband component, which is however essentially irrel- 
evant for the isotropic case in the sense of Anderson's 
theorem. However, intraband scattering by disorder will 
average any anisotropy of the order parameter present in 
the conventional s-wave caso^. In the case of the pnic- 
tides, this intraband component will do the same, with 
the effect of lifting the weak nodes found in our micro- 
scopic calculations. We demonstrate this effect, and its 
consequences, within a simple model where we take weak 
pointlike scatterers, treated in the Born approximation, 
as a model for the intraband part of the true disorder 
scattering in the pnictides. We believe that such scat- 
tering arises primarily from out-of plane dopants like K 
in the hole-doped and F in the electron-doped materi- 
als; since these dopants sit away from the FeAs plane, 
they will certainly produce a significant small-q compo- 
nent of the scattering potential, which will mix states on 
the same sheet. Our results suggest that an extended-s 
wave state with nodes, lifted by disorder in the case of 
the doped superconducting ferropnictides, may explain 
the apparent discrepancies among various measurements 
in the superconducting states of these materials. 
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FIG. 1: (Color online) Fermi surface and extended-s gap 
eigenfunction from Ref. [28l . Fig. 15(d). Red/blue color indi- 
cates sign of order parameter. 



scattering in the superconducting state is 



G(k,a;) 



wTo -I- ekT3 + AkTi 



(2) 



where Cj = lo — Eqi Ck = fk + S3, Ak = Ak + Si, and the 
Sq are the components of the self-energy proportional 
to the Pauli matrices Tq, in particle-hole (Nambu) space. 
If we assume weak scattering, we may approximate the 
self-energy in the Born approximation as 



II. MODEL 



E = n7^|[/(kk')|V3G"(k',c^)T3, 



(3) 



We begin by assuming a superconductor with a sepa- 
rable pair interaction y(k, k') = Vi<I>i((/))<I>i(0'), where 
(j) is an angle parameterizing the electronic momentum k 
on a single circular Fermi surface. To model a situation 
corresponding to an extended-s state on the (3 sheet of 
the Fe-pnictide materials, as found in Ref. H^, we choose 



<I>i((?!)) = 1 + r cos20. 



(1) 



^ 1, ensuring that the pure order 
Ao$i(0) will have nodes near (f) = 



with r = V'/Vi 
parameter Ak = 
7r/4, 37r/4, .... Note that this function is also suitable 
for modelling the d-wave state on the a Fermi surface 
sheets, in the limit r ^ 1. Because the effect of disorder 
will renormalize the constant and cos 2(f> parts of the or- 
der parameter differently, we write it in the general case 
as Ak = Ao -I- A' cos 20. 

The full matrix Green's function in the presence of 



where ni is the concentration of impurities. The self- 
energy has Nambu components 



and 
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, (4) 
, (5) 

• (6) 



As discussed above, we assume pointlike impurity scat- 
tering [/(k, k') = [/qj and treat the disorder in Born ap- 
proximation. We further assume particle-hole symmetry 
such that E3 = 0, leading to Nambu self-energy compo- 
nents after integration perpendicular to the Fermi surface 
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So((/),u;) = r 



Ei((/),a;) 




(7) 



(8) 



where ()0 indicates averaging over the circular Fermi sur- 
face, and r = -kuiNqUq is the normal state scattering 
rate, with A'o the density of states at the Fermi level. 



III. RESULTS 

The BCS gap equation 

Ak = iTrr^^Fk,k'TiG(k',w„) (9) 




then reduces to 



FIG. 2: (Color online) Suppression of the critical temperature 
Tc/Tco vs. normal state scattering rate F/Tco for anisotropy 
parameter r = l.l{dashed), 1.3{solid). The inset shows the 
two curves at smaller values of the scattering rate. The Tc 
suppression curve expected from Abrikosov-Gor'kov theory^ 
(dotted) is included for comparison. 



Ao = iVo [ViAoh + {V'Ao + ViA')l2 + V'A'h] 
A' = No[V'AoIi + {V'A' + rV'Ao)l2+rV'A'l3] 

(10) 

with 

'(cos2(/))™-i\ 




(11) 



and where ujn — {2n + l)7rr is a fermionic Matsubara 
frequency. 



A. Tc suppression 



Near h ^ C , h ^ and I3 = C/2 , where C 



log 



, leading to a critical temperature of 



Tc = l.lSujc cxp 



NoiV,^ + V'y2) 



(12) 



in the clean limit. In the dirty case. Tc is suppressed by 
ordinary disorder until the gap anisotropy is completely 
washed out, at values of the normal state scattering rates 
F many times larger than T^, as also found by Markowitz 
and Kadanofpifl. In Fig. [21 we plot this behavior by solv- 
ing pI7)) with $1 given by Eq. The marginal case 
r = 1 describes the situation where the nodes just touch 
the Fermi surface without disorder, which has been dis- 
cussed earlier—. Note that while the scale of the plot im- 
plies that the critical temperature is rapidly suppressed, 
the scattering rate scale over which this occurs is actu- 
ally much greater than one would expect in. e.g. a d-wave 
superconductor, where a critical concentration Tc Tco 
suffices to suppress superconductivity completely. 
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FIG. 3: (Color online) Normalized spectral gap ^lai^P) /TcO 
vs. angle (j) on the Fermi surface for an extended s-wave state 
with r = 1.3 and F/T^o = (dashed), 0.3 (dotted), 1.0 (solid), 
and 3.1 (dashed-dotted). 



B. Density of states and spectral gap 

In the symmetry broken state, the order parame- 
ter Ak is renormalized by the off-diagonal self-energy 
I]i(k, cj), but the true spectral gap is determined by 
both Nambu components Ei and So. One may define an 
angle-dependent spectral gap flG{4>) by examining the 1- 
particle spectral function A(k, a;) = — (l/7r)ImGii(k, w), 
and plotting either the peak or the energy between the 
peak and the Fermi level where A falls to one-half its 
peak value. We adopt the latter definition here since it 
is similar to one traditionally used in the ARPES com- 
munity, but other measures (e.g. plotting the peak in 
y4(k, w) or simply A(a; = 0,0)) produce very similar re- 
sults. In Fig. [31 we show the s-wave spectral gap plotted 
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FIG. 4: (Color online) Normalized density of states N(uj)/No 
vs. energy uj for the same parameters and line types as Fig. 

El 
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FIG. 5: (Color online) Superfluid density ps/m* vs. T /Tco 
for r — 1.3 and same disorder parameters and line types as in 

Figs, m 



over the Fermi surface for various values of the scattering 
rate F. It is clear that the effect of scattering is as de- 
scribed above: the order parameter is averaged over the 
Fermi surface, which eventually has the effect of lifting 
the nodes, and leading to a true gap in the system, which 
becomes isotropic at sufficiently large scattering rates. 

This can be observed as well in the total density of 
states, which we exhibit in Fig. 21 In the clean situa- 
tion, there are two coherence-type peaks, corresponding 
to the large and small antinodal order parameter scales 
observed in Fig. [31 The addition of disorder smears these 
peaks initially and suppresses the maximum gap feature. 
This is consistent with the expected behavior of a dirty 
nodal superconductor. This behavior continues until the 
point when the nodes are actually lifted, and the system 
acquires a true gap with sharp coherence peaks. When 
the gap has become completely isotropic due to disorder 
averaging, the DOS is identical to the usual BCS DOS, 
as required by Anderson's theorem. 



C. Superfluid density 

The temperature dependence of the superfluid den- 
sity ps (T) reflects the distribution of quasiparticle states 
which contribute to the normal fluid fraction. Early 
data on powdered samples of R-FeAsO^-Fi-a; (R =Pr— , 
Smi ,Nd^) near optimal doping showed exponential T de- 
pendence, as did measurements on Baa;Ki_a;Fe2As2 (Ba- 
122)^. More recent experiments on the latter system 
doped with Co found a power law temperature depen- 
dence close to which evolved towards exponential with 
increasing Co concentration^. A dependence is char- 
acteristic of a dirty system with linear nodes^, but the 
authors of Ref. H concluded the was more likely due 
to nonlocal electrodynamics^, or pairbreaking due to in- 
homogeneity or inelastic scattering. Very recently, how- 
ever, a linear T dependence was measured in the original 
ferropnictide superconductor LaFeFO^^ with Tc — 6K. 



Because this material is stoichiometric, crystals have very 
long mean free paths of more than a thousand Angstrom 
and are capable of supporting dHVA oscillations^. It 
seems very unlikely that anything but order parameter 
nodes can lead to such a power law. The pure s_|_/_ state 
we consider here with weak nodes will immediately yield 
Ps(0) — Ps{T) ~ T, so we take the parameters of the 
previous section and calculate ps directly. 

Within the current BCS-type model, the xx compo- 
nent of the superfluid density tensor ps may be written 



cos 



duj tanh — Re ■ 
2T 



- A2 )3/2 



(13) 

where p is the full electron density. Note this stan- 
dard expression for the superfluid density of a dirty 
superconductor— is complete within mean field theory 
since vertex corrections to the current response vanish 
for singlet superconductors and pointlike scatterers^. In 
Fig. IHl we show the superfluid density = (ps.xx + 
Ps.yy)/2 calculated from Eq. (|13p vs. temperature for 
various scattering rates^. In the pure system, the su- 
perfluid density is linear at the lowest temperatures, re- 
flecting the existence of line nodes. At a temperature 
corresponding to the smaller antinodal gap^, there is a 
decrease in the rate at which thermally excited quasipar- 
ticles depopulate the condensate, and therefore a change 
in slope as visible in the figure. We note that this kink- 
like behavior may disappear when excitations from other 
Fermi surface sheets which contribute at higher temper- 
atures are included, but also that this type of upward 
curvature in ps vs. T was observed by Fletcher et ali^ in 
the cleanest LaFePO samples. 

If disorder is now added without lifting the nodes, the 
generic behavior will be quadratic in temperature^, i.e. 
Ps — /5s(0)(l — aT^/F^), where a is a constant of order 
unity. When the nodes are completely lifted, an expo- 
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FIG. 6: (Color online) logio(l - p,(r)/p.(0)) vs. log^j T/Tc 
for r — 1.3 and various scattering rates (Circles: F = 0, 
triangles: T/Tco = 0.3, crosses: F/TcO = F, and squares: 
r/Tco = 3.1 . Solid line: T; dashed-dotted line: T^. 



ncntial T dependence, ps — Ps(0)(l — acxp(— A,„i„/T)) 
must dominate at the lowest temperatures. These power 
laws are displayed more precisely in the log-log plot 
of Fig. [51 The pure system follows a linear-T law 
(ps(0) — Ps{T) ~ T) down to the lowest temperatures, as 
expected because of the linear nodes, whereas the dirty 
systems follow a law over an intermediate T range, 
also expected because the states are uniformly broad- 
ened in this range by Born scattering with 1/r roughly 
linear in energy. At lower temperatures, however, the 
dirty cases no longer exhibit power laws in temperature, 
with the curves showing fully developed gaps in Fig. 3] 
following activated behavior. The T/Tco = 0.3 case is 
marginal in the sense that the effective spectral gap has 
nodes which just touch the Fermi surface, and the figure 
shows that when the system is close to this condition the 
low-T behavior is not a simple power law or exponential. 



IV. EFFECT OF ADDITIONAL FERMI 
SURFACE SHEETS 

While the desired effect of node lifting by disorder 
for an extended-s wave state with weak nodes has now 
been exhibited, the amount of disorder required to lift 
the nodes has also been shown to suppress Tc substan- 
tially. As strong sample-to-sample variations of Tc have 
not been observed in the ferropnictidcs, this seems in- 
consistent with the current overall body of experimental 
data. There are several possibilities to explain this dis- 
crepancy. The first is that the nodes, e.g. in the LaFePO 
system, may be accidentally extremely weak. This seems 
unlikely, since the linear-T behavior extends over a sig- 
nificant fraction of Tc in this material. Were the nodes 
to be marginal or extremely weak, the phase space from 
near-quadratic nodes would actually lead to a T^/^ tem- 
perature variation of the penetration depth in the pure 



FIG. 7: (Color online) Spectral gaps QG{(f))/Tco vs. 4> for the 
a (top) and (3 (bottom) sheets for the 2-band model. Here 
the pairing parameters Vii = 1.0, V12 ~ —0.6, V22 =1.5 were 
chosen so that r^ff = 1.3 is the same as in the 1-band case. 
The densities of states ratio was taken as A''2(0)/A^i(0) — 1.25. 
Scattering rates for intraband scattering correspond to F = 
7ii7riVi(0)|(7iii^ = (dashed), 0.3 (dashed-dotted), 0.6 (solid). 

system, which is not observed. We have also examined 
the above model with smaller values of r ^ 1, and find 
that substantial scattering rates are still required to cre- 
ate a significant spectral gap. 

It is possible that our restriction to Born (weak) scat- 
tering is inadequate. Unitary scatterers modify primarily 
the states near the nodes, creating a gap without affect- 
ing the states at higher energies; smaller normal state 
scattering rates can therefore produce a large node lift- 
ing effect without suppressing Tc significantly^. It may 
be interesting to explore this effect further, but a priori 
there is no obvious reason for the out-of-plane impurities 
which appear to dominate the doped pnictide materials 
to produce such strong scattering potentials. 

A final possibility, which we do pursue here, is that in 
the ferropnictidcs, the single Fermi surface sheet we have 
studied thus far is coupled to other sheets which con- 
trol the Tc suppression. This appears plausible within 
the context of our microscopic spin-fluctuation pair- 
ing calculations^^, which show that for extended-s type 
states, while the /3 sheets appear to have weak nodes, 
the a sheets are more isotropic for the extended-s wave 
solutions. The pairing on the a and /3 sheets is strongly 
coupled by the tt, (unfolded zone) scattering processes 
which dominate the spin-fluctuation spectrum. To sup- 
press the critical temperature, it is clear that not only 
the (3 condensate, but also the a condensate must be sup- 
pressed. The a condensate may therefore be expected to 
act as a reservoir which maintains the critical tempera- 
ture (provided interband scattering by disorder is rela- 
tively weak), while the /3 condensate is nodal and there- 
fore dominates low-temperature and low-energy proper- 
ties. 

To test this hypothesis within the philosophy of the 
current paper, we take the above model and add to it a 
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single additional sheet with constant pairing amplitude, 
and couple the two sheets by constant pairing potentials, 
leading to a total pairing interaction. 



V{k,k') = Fi$i(k)$i(k')+t^2$2(k)$2(k') 

+Vi2[M^)M^') + <^2(k)0i(k')], (14) 

where (pi (k) = 1 + r cos 20 is understood as before to 
describe states with k on the /3 sheet, whereas (^2(k) = 1 
for k on the a sheet. Note that V12 is chosen of oppo- 
site sign to Vi,V2 so as to induce a sign-changing order 
parameter between the two sheets. 

We now consider the scattering from nonmagnetic im- 
purities. As discussed by Refs. [33''34','35', it is convenient 
to parameterize the scattering potential in terms of an 



amplitude U22 describing scattering k 



the 



sheet, similarly Un on the /3 sheet, and t/12,21 between 
the two, such that the Born self-energy becomes 



E(kG l,uj) 



ni 



k'ei 



J2\U,2\'Gik\u;) 



k'e2 



(15) 



and similarly for k G 2. To compare the 1-band and 
2-band cases, we now choose pairing parameters such 
that the order parameter on the (3 sheet is nearly iden- 
tical to that we had in the 1-band case above with 
reff = (Amax " Amin) / {^Tnax + Amin) = 1-3. This is il- 
lustrated as the pure spectral gap ilc in Fig- where the 
corresponding isotropic gap on the a sheet is also shown. 
As disorder is added, the nodes are removed from the /3 
sheet, as before, and a full gap created. At the same time 
the a gap is only very slightly suppressed in this process. 
These features are evident also in the total DOS shown 
in Fig. [5] for the same scattering parameters. 

Our aim in this section is to see if within a two-band 
picture one can understand qualitatively under what cir- 
cumstances a substantial spectral gap can be opened on 
the Fermi surface without suppressing Tc substantially 
in the simple 1-band example. The naive hypothesis for- 
mulated above, that the presence of a more isotropic gap 
on a second Fermi surface sheet is sufficient to "protect" 
Tc against intraband scattering, is not universally cor- 
rect. This is because the order parameters on the two 
sheets are strongly coupled by the interband pair inter- 
action. To illustrate which aspects of the problem are 
important for the relative proportion of gap creation rel- 
ative to Tc suppression, we compare in Fig. [9] the gap 
created by a certain amount of disorder versus the cor- 
responding suppression of Tc, for the 1-band case and 
two 2-band cases where the effective 1 sheet anisotropy 
reff is held fixed, but the ratio of the densities of states 
on the two sheets are varied. It is seen that the critical 
temperature is rendered more robust when the density of 
states A''2(0), which controls the pairing weight on sheet 2 




FIG. 8: (Color online) Density of states N{u) / Nxotai vs. 
oj/Tco for 2-band system with re// = 1.3. Same parame- 
ters and line types as in Fig. [Tjlnsert: expanded low-energy 
region. 
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FIG. 9: (Color online) Normalized spectral gap 
m\n^G{(t>) /Tco vs. Tc suppression 1 — Tc/Tao for vari- 
ous impurity concentrations comparing 1-band and 2-band 
models. Parameters are chosen such that re// = 1.3 on 
the 1 (/3) sheet for the pure superconducting state in all 
cases. Diamonds: 1-band model with Tc/Tco and Q.g/Tco 
taken from Figs. [2] and [S] respectively. Squares: two-band 
model with Vii = 1.0, V12 = —0.6,1/22 =1.5, with densities 
of states ratio A^2(0)/A^i(0) = 1.25 and 1 sheet anisotropy 
parameter r = 1.76. Circles: two-band model with same 
Vii = 1.0, V12 = —0.6, V22 =1.5, but with densities of states 
ratio iV2(0)/iVi(0) = 1.5 and r = 1.88. 



(a), is increased. Any effect which enhances the nodeless 
sheet 2 pairing weight makes Tc less susceptible to intra- 
band scattering. Thus we conclude that the amount of 
Tc suppression associated with gap creation depends on 
the details of the situation, and can become quite small. 

The true Fermi surface structure is of course more com- 
plicated than the 2-sheet model considered here, so it is 
interesting to ask whether results from the more com- 
plete model show the desired effect. We show in Fig. 10 
calculations for the spectral gap in the extended s-wave 
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FIG. 11: (Color online) minSlG/TcO vs. 1 -TrJTco for 
disorder parameters as Fig. 1101 



same 
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FIG. 10: (Color online) Spectral gap vs. Fermi surface arc 
length in arbitrary units on each of four Fermi sheets (ai, a2, 
/3i arid B2) for realistic spin fluctuation model of Fig. 17 of 
Ref. [2^. Curves correspond to an arbitrary range of impu- 
rity concentrations, from red (clean) to blue (dirty). Bottom 
panel: detail of nodal region of /3 sheets. 



{Aig) state obtained from the microscopic calculations 
of Rcf. m for various values of intraband disorder scat- 
tering. The node lifting phenomenon on the (3 sheets is 
clearly observed. In Fig. [Til we show that relatively little 
Tc suppression accompanies this realistic case; significant 
gaps of order 0.1 — Q.2Tc are obtained for 1 — Tc/Tca of 
only ~ 10%. 



V. CONCLUSIONS 

There is now evidence that the superfluid density has 
a T dependence consistent with a fully developed gap in 
some samples, while power laws in T, including linear 
T, have been reported in others. It is possible that pa- 
rameters related to the electronic structure of the pure 
state tune the various materials such that different su- 
perconducting ground states are realized. This is "natu- 



ral" due to the proximity of these systems to a situation 
where nearly circular Fermi surface sheets nest perfectly, 
in which case spin fluctuation theory predicts a degener- 
acy between the extended-s {Aig) and d {Big) states. 

Here we have explored an alternative possibility which 
assumes that there are nodes in an cxtended-s-wave 
{Aig) state for the ideal ferropnictide material, but that 
these nodes are lifted by small momentum disorder scat- 
tering. In this regard, it is interesting that the stoichio- 
metric ferropnictide LaFePO (with a mean free path of 
more than lO^nm and capable of supporting dHvA oscil- 
lations) displays a linear T dependence of the superfluid 
density at low temperatures. While many mechanisms 
can give rise to a power law T dependence in the super- 
fluid density, we are aware of only one explanation for a 
linear-T law, namely the existence of nodes in the pure 
state. On this basis we speculate that the exponential 
behaviors observed in other materials are due to disor- 
der. On the other hand, evidence for the opposite trend 
has been reported in the Bai_2;K:i.Fe2As2 system, where 
the cleanest crystals appear to exhibit exponential T de- 
pendence, whereas more disordered samples exhibit a low 
temperature dependence which mimics a power law. 

We have therefore explored a scenario in which weak 
nodes in a sign-changing s_|_/_ state on the (3 sheets of 
the ferropnictides are "lifted" by nonmagnetic disorder. 
Qualitative aspects of this phenomenon are obtainable 
already within a simple 1-band model where weak point- 
like scatterers are accounted for. With increasing dis- 
order, the nodes disappear and are replaced by a fully 
developed gap and an activated temperature dependence 
of ps. Thus one could imagine that the entire class of 
ferropnictide materials has intrinsic nodes in the ground 
state order parameter of the analog pure system, which 
then disappear in most of the doped, much dirtier ma- 
terials. We have shown that in the simplest model the 
lifting of the nodes corresponds to a significantly larger 
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pairbreaking rate and concomitant suppression than 
observed experimentally, but that this undesirable aspect 
of the theory is substantially eliminated by the addition 
of additional bands with more isotropic pair state, as 
found in the microscopic theory2^. 

The justification for our neglect of intcrband scattering 
in most of this work is based on the empirical robustness 
of Tc to differences in sample quality, suggesting that 
intcrband scattering in the sign-changing s state is neg- 
ligible. In addition, Yukawa-type models of realistic im- 
purity potentials with finite range have small large-g am- 
plitudes. Nevertheless, for a quantitative description of 
these materials both types of scattering should be consid- 
ered, and a more microscopic understanding of the types 
of potentials introduced, e.g. by a K replacing Ba,Sr out 
of the FcAs plane and by a Co substituting directly for 
an Fc would be very useful. 

The simple theory presented in Section [IT] is actually 
applicable to a state with arbitrary anisotropy parame- 
ter r, allowing us to describe the pure d-wave state for 
r ^ 1, as well as a state with deep minima of the order 
parameter ( "quasi-nodes" ) , but no true nodes for r < 1 
as found, e.g. in Ref. [2^. Thus we emphasize that it is 
possible that small differences in the electronic structure 
of different materials may give rise to slightly different 
gap structures, and even different symmetry classes^^ in 
different materials. The implication of current ARPES 
experiments, that the order parameter anisotropy around 
the various Fermi surface sheets is absent or very weak 
is, however, difficult to reconcile with experimental data 
which indicate low energy excitations, and a disorder- 



based explanation of the kind given here, where nodes or 
quasi-nodes are lifted by momentum averaging, seems to 
us the most likely way to understand the existing data. 
This picture will have implications for other properties, 
including transport properties, as discussed in Ref. |5l| 
for the 1-band situation, but there may be interesting 
modifications peculiar to the multiband case. Work on 
this direction is in progress. 

We note in closing that an alternative approach takes 
the point of view that the intrinsic material is fully 
gapped with a sign-changing s-wave state, and that the 
intcr-Fcrmi surface scattering leads to the appearance of 
low-temperature power law behavio r'^^i^^i'^^ . Clearly sys- 
tematic electron irradiation or some other source of dis- 
order which does not change the doping of the system 
would be a very important experimental way to distin- 
guish this proposal from our own. 
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